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A 3-fluid hydrodynamic model is introduced for simulating heavy-ion collisions 
I at incident energies between few and about 200 A-GeV. In addition to the two 
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baryon-rich fluids of 2-fluid models, a new model incorporates a third, baryon-free 
(i.e. with zero net baryonic charge) fluid which is created in the mid-rapidity region. 



CN ' fluid neither thermalizes nor interacts with the baryon-rich fluids. After formation it 

0^ , thermalizes and starts to interact with the baryon-rich fluids. It is found that for r 



the interaction strongly affects the baryon-free fluid. However, at reasonable finite 
formation time, r ~ 1 fm/c, the effect of this interaction turns out to be substantially 
reduced although still noticeable. Baryonic observables are only slightly affected by 
the interaction with the baryon-free fluid. 



I. INTRODUCTION 

Nearly fifty years have passed since the paper "Relativistic kinetic equation" by Spar- 
tak T. Belyaev and Gersh I. Budker was published The proposed relativistic formula- 
tion of the distribution function and the kinetic equation with small-angle scattering have 
been included in many textbooks and found numerous applications in atomic physics and 
electron-positron plasma. Recently a generalized relativistic kinetic equation of this type 
was implemented for describing the partonic evolution in very early stages of a heavy-ion 
collision at ultra-relativistic RHIC energies j^. In the present paper we address more mod- 
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erate, but nevertheless highly relativistic, energies, i.e. up to those reached at the CERN 
SPS. The relativistic kinetic equation is used in a peculiar way, namely to derive a coupling 
term for three-fluid hydrodynamic equations. 

Two-fluid hydrodynamics with free-streaming radiation of pions was advanced first in 0, 
The initial stage of heavy-ion collisions definitely is a highly non-equilibrium process. 
Within the hydrodynamic approach this non-equilibrium is simulated by means of a 2-fiuid 
approximation, which takes care of the finite stopping power of nuclear matter 0,0], and 
simultaneously describes the entropy generation at the initial stage. The radiated pions form 
a baryon-free matter in the mid-rapidity region, while two baryon-rich fluids simulate the 
propagation of leading particles. The pions are the most abundant species of the baryon- 
free matter which may contain any hadronic and/or quark-gluon species including baryon- 
antibaryon pairs. 

First applications of the 2-fluid model [3, to the description of heavy-ion collisions in the 
wide range of incident energies, from those of SIS to SPS, were quite successful. One of the 
advantages of the hydrodynamic models is that they directly address the equation of state 
(EoS) of nuclear matter, which is of prime interest for this domain of physics. In particular, 
we have shown recently that the experimental excitation function of the directed flow is 
well described by the mixed-phase EoS in contrast to earlier predictions of the two-phase 
bag-model EoS. In these 3D hydrodynamic simulations we describe the whole process of 
the reaction, i.e. the evolution from the formation of a hot and dense nuclear system to its 
subsequent decay. This is in distinction to numerous other simulations, which treat only 
the expansion stage of a fireball formed in the course of the reaction, while the initial state 
of this dense and hot nuclear system is constructed from either kinetic simulations or more 



general albeit model-dependent assumptions (e.g. see |lOl lllj). 

However, the approximation of free-streaming pions, produced in the mid rapidity region, 
was still irritating from the theoretical point of view, in particular, because the relative mo- 
menta of the produced pions and the leading baryons are in the range of the A resonance 
for the incident energies considered. This would imply that the interaction between the 
produced pions and the baryon-rich fiuids should be strong. The free-streaming assumption 
relies on a long formation time of produced pions. Indeed, the proper time for the formation 
of the produced particles is commonly assumed to be of the order 1 fm/c in the comov- 
ing frame. Since the main part of the produced pions is quite relativistic at high incident 
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energies, their formation time should be long enough in the reference frame of calculation 
to prevent them from interacting with the baryon-rich fluids. However, this argument is 
qualitative rather than quantitative, and hence requires further verification. The first at- 
tempt to do this was undertaken by the Frankfurt group which started to explore an 
opposite extreme. They assumed that the produced pions immediately thermalize, forming 
a baryon-free fluid (or a "fireball" fluid, in terms of and interact with the baryon-rich 
fluids. No formation time was allowed, and the strength of the corresponding interaction 
was guessed rather than microscopically estimated. This opposite extreme, referred to as a 
(2-|-l)-fluid model and being not quite justified either, yielded results substantially different 
from those of the free-streaming approximation. This was one of the reasons why in subse- 
quent applications the Frankfurt group neglected the interaction between baryon-free and 
baryon-rich fluids while keeping the produced pions thermalized thus effectively restor- 
ing the free-streaming approximation. However, the assumed immediate thermalization of 
the flreball fluid together with the lack of interaction with baryon-rich fluids still was not a 
consistent approximation. 

In this paper we would like to return to the problem of veriflcation of the free-streaming 
approximation for the produced pions. To do this, we extend the 2-fluid model of refs. 

to a 3-fluid model, where the created baryon-free fluid (which we call a 
"flreball" fluid, according to the Frankfurt group) is treated on equal footing with the baryon- 
rich ones. However, we allow a certain formation time for the flreball fluid, during which 
the constituents of the fluid propagate without interactions. Furthermore, we estimate the 
interaction between flreball and baryon-rich fluids by means of relativistic kinetic equation 
and elementary cross sections. 

In this paper we consider incident energies in the range range from few to about 200 
A-GeV (i.e. from AGS to SPS energies). The interest to this this energy was recently revived 
in connection with the project of the new accelerator facility at GSI SIS200 The goal of 
the research program on nucleus-nucleus collisions at this planned facility is the investigation 
of nuclear matter in the region of incident energies (-Eiab — 10 — 40 A-GeV), in which 
the highest baryon densities and highest relative strangeness at moderate temperatures are 
expected. Here, the QCD phase diagram is much less explored, both experimentally and 
theoretically, as compared to the higher energy region characterized by higher temperatures, 
but lower net baryon densities, where lattice QCD calculations [l^ and a large body of 
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experimental data from SPS (CERN) and RHIC (BNL) are available by now. 



II. 3-FLUID HYDRODYNAMIC MODEL 

Unlike the one-fluid hydrodynamic model, where local instantaneous stopping of projec- 
tile and target matter is assumed, a specific feature of the dynamic 3-fiuid description is a 
finite stopping power resulting in a counter-streaming regime of leading baryon-rich mat- 
ter. Experimental rapidity distributions in nucleus-nucleus collisions support this counter- 
streaming behavior, which can be observed for incident energies between few and 200 A-GeV. 
The basic idea of a 3-fiuid approximation to heavy-ion collisions 0, 0, is that at each 
space-time point x = (t, x) the distribution function of baryon-rich matter, /br(a;,p), can be 
represented as a sum of two distinct contributions 

fhr{x,p) = fp{x,p) + /t(x,p), (1) 

initially associated with constituent nucleons of the projectile (p) and target (t) nuclei. 
In addition, newly produced particles, populating the mid-rapidity region, are associated 
with a fireball (f) fluid described by the distribution function ff{x,p). Note that both the 
baryon-rich and fireball fluids may consist of any type of hadrons and/or partons (quarks 
and gluons), rather then only nucleons and pions. However, here and below we suppress the 
species label at the distribution functions for the sake of transparency of the equations. 

With the above-introduced distribution functions /„ (a =p, t, f), the coupled set of 
relativistic Boltzmann equations looks as follows: 

p,d^U{x,p) = Cp(/p,/0 + Cp(/p,/f), (2) 
p,d^Mx,p) = a(/p, /t) + a(/t, /f), (3) 

p,d^ji{x,p) = a(/p, A) + Q(/p, /f) + /f), (4) 

where Ca denote collision terms between the constituents of the three fluids. We have 
omitted intra-fluid collision terms, like Cp(/p,/p), since below they will be canceled any 
way. The displayed inter-fluid collision terms have a clear physical meaning: Cp/j(/p,/t), 
Cp/t{fp/t, ff), and Ci{fp/t,fi) give rise to friction between p-, t- and f-fluids, and Cf(/p, /t) 
takes care of particle production in the mid-rapidity region. Note that up to now we have 
done no approximation, except for hiding intra-fluid collision terms. 
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Let us proceed to approximations which justify the term "fluids" having been used al- 
ready. We assume that constituents within each fluid are locally equilibrated, both thermo- 
dynamically and chemically. In particular, this implies that the intra-fluid collision terms 
are indeed zero. This assumption relies on the fact that intra-fluid collisions are much more 
efficient in driving a system to equilibrium than inter-fluid interactions. As applied to the 
flreball fluid, this assumption requires some additional comments, related to the concept of 
a flnite formation time. During the formation proper time r after production, the flreball 
fluid propagates freely, interacting neither with itself nor with the baryon-rich fluids. . Af- 
ter this time period, the flreball matter locally thermalizes and starts to interact with both 
itself and the baryon-rich fluids. Being heated up, these three fluids may contain not only 
hadronic and but also deconflned quark-gluon species, depending on the EoS used. 

The above assumption suggests that interaction between different fluids should be treated 
dynamically. To obtain the required dynamic equations, we flrst integrate the kinetic 
Eqs. over momentum and sum over particle species with weight of baryon charge. 

This way we arrive to equations of the baryon charge conservation 

d,Ji:{x) = 0, (5) 

for a =p and t, where Jj^ = UaU'^ is the baryon current deflned in terms of baryon density 
ria and hydrodynamic 4-velocity normalized 1. Eq. © implies that there is 

no baryon-charge exchange between p- and t-fluids, as well as that the baryon current of 
the flreball fluid is identically zero, = 0. Integrating kinetic Eqs. over momentum 

with weight of 4-momentum p'^ and summing over all particle species, we arrive at equations 
of the energy-momentum exchange for energy-momentum tensors T^'^ of the fluids 

W^(^) = -F;ix)+K(^), (6) 

d.rrix) = -Fax) + Fli{x), (7) 
d.rrix) = F;{x) + Fax) - F^^{x) - F^xl (8) 

where the F^ are friction forces originating from inter-fluid collision terms in the kinetic 
Eqs. ((21)^0 • F^ and F^ in Eqs. (jHI)^© describe energy-momentum loss of baryon-rich 
fluids due to their mutual friction. A part of this loss \F^ — F^ is transformed into thermal 
excitation of these fluids, while another part (Fp -|- F^) gives rise to particle production into 
the flreball fluid (see Eq. (jH}). F^ and F^^^ are associated with friction of the flreball fluid with 
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the p- and t-fluids, respectively. Note that Eqs. ©^(0) satisfy the total energy-momentum 
conservation 



0. 



(9) 



As described above, the energy-momentum tensors of the baryon-rich fluids {a =p and 
t) take the conventional hydrodynamic form 



rpfJLV 

a 



(e„ + Pa)< K 



(10) 



in terms of the proper energy density, Ea, and pressure, P^. For the fireball, however, only 
the thermalized part of the energy-momentum tensor is described by this hydrodynamic 
form 



(11) 



Its evolution is defined by a Euler equation with a retarded source term 



^^yK)M^'(^) = / d^x'S' [x-x'- Uf{x)t) 



F;{x')+Fax' 



K{x)-Fr,{x), 



[12) 



where r is the formation time, and 



F;{x') + Fax') 
\FJx')+FAx') 



(13) 



is a free-streaming 4- velocity of the produced fireball matter. In fact, this is the velocity 
at the moment of production of the fireball matter. According to Eq. ()12|). the energy 
and momentum of this matter appear as a source in the Euler equation only later, at the 
time UpT after production, and in different space point x' — \Jf{ compared to the 

production point x'. From the first glance, one can immediately simplify the r.h.s. of Eq. 
(fT^ by performing integration with the 5-function. However, this integration is not that 
straightforward, since the expression under the ^-function, x — x' — Uf{x')t = 0, may have 
more than one solution with respect to x'. The latter would mean that the matter produced 
in several different space-time points x' simultaneously thermalizes in the same space-time 
point X. This is possible due to the nonlinearity of the hydrodynamic equations. 
The residual part of T/"^ (the free-streaming one) is defined as 



(cq) fiu 



(14) 
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The equation for T^^'^^^'^ can be easily obtained by taking the difference between Eqs. (jH)) and 
()12|). If all the fireball matter turns out to be formed before freeze-out, than this equation 
is not needed. Thus, the 3-fluid model introduced here contains both the original 2-fluid 
model with pion radiation 1^,0,0,0,0 and the (2+l)-fiuid model Q, 13] as limiting cases 
for r — > oo and r = 0, respectively. 

The nucleon-nucleon cross sections at high energies are strongly forward-backward 
peaked. In this case the Boltzmann collision term can be essentially simplified, since the in- 
volved 4-momentum transfer is small. The small-angle-scattering expansion of the collision 
integral results in the relativistic Fokker-Planck equation, as first derived by Belyaev and 

n n 

Budker Precisely this equation was used in p to estimate the friction forces, and 
F]\ proceeding from only A^A^ elastic scattering. Later these friction forces were calculated 
19| based on (both elastic and inelastic) experimental inclusive A^A^ cross sections 

F: = PpPt [iK - ur^ Dp + [ul + <) De] , (15) 

a =p and t, p =t and t =p. Here, pa denotes the scalar densities of the p- and t-fluids, 

Dp/E = friN Kci (^P/E{spt), (16) 

where is the nucleon mass, Spt = (up + Wt) is the mean invariant energy squared of 
two colliding nucleons from the p- and t-fluids, V^^\ = [spt(spt — 4m^)]^/^/2m|^ is the mean 
relative velocity of the p- and t-fluids, and <Jp/E{spt) are determined in terms of nucleon- 
nucleon cross sections integrated with certain weights (see jsi J, Q, 0,13 for details). It was 
found in 19] that a part of these friction terms, which is related to the transport cross-section, 
may be well parameterized by an effective deceleration length Aefr with a constant value 
Aeflf ~ 5 fm. However, there are reasons to consider Aeg as a phenomenological parameter, 
as it was pointed out in 7]. Indeed, as it is seen from Eq. (jlfjj) . this friction is estimated 
only in terms of nucleon-nucleon cross sections while the excited matter of baryon-rich 
fluids certainly consists of great number of hadrons and/or deconfined quarks and gluons. 
Furthermore, these quantities may be modified by in-medium effects. In this respect, -Dp/E 
should be understood as quantities that give a scale of this interaction. 

Eqs. (jni)^© and (fT^ . supplemented by a certain EoS and expressions for friction forces 
F", form a full set of equations of the relativistic 3-fluid hydrodynamic model. To make this 
set closed, we still need to define the friction of the fireball fluid with the p- and t-fluids, F^ 
and F^^ in terms of hydrodynamic quantities and some cross sections. 
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III. INTERACTION BETWEEN FIREBALL AND BARYON-RICH FLUIDS 

Our aim here is to estimate the scale of the friction force between the fireball and baryon- 

I I 

rich fluids, similar to that done before for baryon-rich fluids [19]. To this end, we consider 
a simplified system, where all baryon-rich fluids consist only of nucleons, as the most the 
abundant component of these fluids, and the fireball fluid contains only pions. 

For incident energies from 10 (AGS) to 200 A-GeV (SPS), the relative nucleon-pion en- 
ergies are in the resonance range dominated by the A-resonance. To estimate this relative 
energy we consider a produced pion, being at rest in the center of mass of the colliding 
nuclei, = {m^, 0, 0, 0}cm- Baryon-rich fluids decelerate each other during their interpen- 
etration. This means that the nucleon momentum should be smaller than the incident 
momentum, \p^\ < \{miy'ycm,Pcm}\, where is the gamma factor of the incident nucleon 
in the cm frame. Calculating the invariant relative energy squared s = {p + q)"^ at -Eiab = 158 
A-GeV, we obtain s^/^ < 1.8 GeV. This range of s precisely covers the resonance region, 1.1 
GeV < 5^/2 < 13 GeV ^20]. At Ei^b = 10 A-GeV we arrive at s^/^ < 13 Qgy, which is also 
within the resonance region. At even lower incident energies the strength of the fireball fluid 
becomes so insignificant, as compared with thermal mesons in the p- and t-fluids, that the 
way of treatment of its interaction with the baryon-rich fluids does not essentially affect the 
observables. For the same reason we do not apply any special prescription for the unification 
of the fireball fluid with the baryon-rich fluids, since this may happen only at relatively low 
incident energies Ei^h < 10 A-GeV. 

The resonance-dominated interaction implies that the essential process is absorption of 
a fireball pion by a p- or t-fluid nucleon with formation of an i?-resonance (most probably 
A). This produced i?-resonance still belongs to the original p- or t-fluid, since its recoil due 
to absorption of a light pion is small. Subsequently this i?-resonance decays into a nucleon 
and a pion already belonging to the original p- or t-fluid. Symbolically, this mechanism can 
be expressed as 

AT" + vr^ ^ i?" ^ AT" + tt". 

As a consequence, only the loss term contributes to the kinetic equation for the fireball fluid. 
Proceeding from the above consideration, we write down the collision term between 
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fireball-fluid pions and a-fluid nucleons (a =p or t) as follows 

a(/.,/f) = - iV^^-^is) ft''\x,p) faix,q), (17) 

where s = (p + g)^, 

is the rate to produce a baryon i?-resonance, and cr^ot^^i^) parameterization of 

experimental pion-nucleon cross-sections j^]. Here, only the distribution function of formed 
(and hence thermalized) fireball pions, f^^'^\ enters the collision term, since the non- formed 
particles did not participate in the interaction by assumption. 

Integrating C{{fa, fi) weighted with the 4- momentum p'^ over momentum, we arrive at 

J Qo Po 

- ^^^0 [J [J —P P = A.^P« , (18) 

where we substituted p° and s by their mean values, < p° >= m.„Uf and Sfa = {fnirUi + 
mjqUaY, and introduced the transport coefficients 

D,^ = l^^"^^(sfa)/(m;vmO = V,':, a^.r^'M- (19) 

Here, V^^" = [(sfa — — m^)^ — 4m^m^]"'^/^/(2mjv'^7r) denotes the mean invariant relative 
velocity between the fireball and the a-fiuids. Thus, we have expressed the friction in 
terms of the fireball-fiuid energy- momentum density T^'^, the scalar density pa of the a fluid, 
and a transport coefficient Dfa- Note that this friction is zero until the fireball pions are 
formed, since Tf^'^^^''' = Q during the formation time r. 

In fact, the above treatment is an estimate of the friction terms rather than their strict 
derivation. This peculiar way of evaluation is motivated by the form of the final result ()18p. 
An advantage of this form is that m,r and any other mass do not appear explicitly, and 
hence allows a natural extension to any content of the fluid, including deconfined quarks 
and gluons, assuming that Dfa represents just a scale of the transport coefficient. 



IV. SIMULATIONS OF NUCLEUS NUCLEUS COLLISIONS 



The relativistic 3D code for the above described 3-fluic 
of modifying the existing 2-fluid 3D code of refs. jsi 14 



model was constructed by means 
In actual calculations 



10 



we used the mixed-phase EoS developed in j2l|, |22, |23| . This phenomenological EoS takes 



into account a possible deconfinement phase transition of nuclear matter. The underlying 
assumption of this EoS is that unbound quarks and gluons may coexist with hadrons in 
the nuclear environment. In accordance with lattice QCD data, the statistical mixed-phase 
model describes the first-order deconfinement phase transition for pure gluon matter and 
crossover for that with quarks 21 , 2^ [2^ . 



We performed simulations of nucleus-nucleus collisions Pb+Pb at i^iab = 158 A-GeV and 
Au+Au at i^iab = 10.5 A-GeV. General dynamics of heavy-ion collisions is illustrated in 
Fig. 1 by the energy-density evolution of the baryon-rich fluids (eb = £p + ^t, in the cm 
frame of colliding nuclei) in the reaction plane of the Pb+Pb collision. Different stages of 
interaction at relativistic energies are clearly seen in this example: Two Lorentz-contracted 
nuclei (note the different scales along the x- and z-axes in Fig.l) start to interpenetrate 
through each other, reach a maximal energy density by the time ~ 1.1 fm/c and then 
expand predominantly in longitudinal direction forming a "sausage-like" freeze-out system. 
At this and lower incident energies the baryon-rich dynamics is not really disturbed by the 
fireball fluid and hence the cases r = and 1 fm/c turned to be indistinguishable in terms 
of ^b- 

In Fig. 2 the dynamic evolution of the fireball energy density {sf, in the cm frame of 
the colliding nuclei) in the reaction plane of the Pb+Pb collision at impact parameter b = 
2 fm is shown for two values of the formation time, r = fm/c (the left column of panels) 
and 1 fm/c (the right column of panels). It starts to form near the time moment, when the 
maximal energy density ^b is reached. The /-fiuid evolution indeed looks like that for an 
expanding fireball, its density essentially depends on the formation time. 

To quantitively reveal the role of the interaction between fireball and baryon-rich fiuids, 
we followed the evolution of the total energy released into the fireball fiuid 

^(released) = dt' J d^'x' [F^{x') + F,\x')) , (20) 

cf. Eq. dS}, and the total energy kept in the fireball fiuid (both thermalized and nonther- 
malized) after interaction. 



Ei'°'\t) = Jd'xT^^t,^) 



dt' / rfV F°(x') + F°(x') - F° (x') - F°(a:') , (21) 
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Figure 1: Time evolution of the energy density, £b = + ^t) for the baryon-rich fluids in the 
reaction plane {xz plane) for the Pb+Pb collision (E'lab = 158 A-GeV) at impact parameter b = 
2 fm. Shades of gray represent different levels of £b as indicated at the right side of each panel. 
Numbers at this palette show the values (in GeV/fm'^) at which the shades change. Arrows 
indicate the hydrodynamic velocities of the fluids. 
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Figure 2: The same as in Fig. 1, but for the firebah-energy density (sf in the cm frame of the 
cohiding nuclei) for two formation times, r = fm/c (the left column of panels) and 1 fm/c (the 
right column of panels). 



cf. Eq. (jS)), in the cm frame of two colliding nuclei. Results of the calculation are pre- 
sented in Fig. 3. To provide a common scale, the quantity calculated with the 
formation time r = 100 fm/ c is presented in all the panels. The r = 100 fm/ c case practically 
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implies absence of interaction between the fireball and baryon-rich fluids and the equality 
^(tot) _ ^(released) ^ bccausc t/^^^^"" = and hcucc = = 0, cf. Eq. (HHD- 




5 10 5 10 6 10 

t [fm/c] t [fm/c] t [fm/c] 



Figure 3: Time evolution of the total energy released into the fireball fluid, j^'^^'^'^^'^ (solid lines), 
and the total energy kept in the fireball fiuid after interaction, E^^°^^ (dashed-dotted lines). Two 
nucleus-nucleus collisions Pb+Pb at Fiab = 158 A-GeV (three upper panels) and Au+Au at i?iab = 
10.5 (three bottom panels) A-GeV, both at zero impact parameter, were calculated with different 
formation times r indicated in the panels. The upper dashed curve in all 6 panels represent 
^(released) (,j^jg^jg^|;gjj with the formation time r = 100 fm/c. 



First, we see that the energy release into in the fireball fluid occurs only during a short 
time of interpenetrating of colliding nuclei. As it has been expected, at zero formation time 
r = the interaction with the baryon-rich fluids strongly affects the fireball fluid: It reduces 
its total energy Fj*"*-* as compared to the case without interaction (i.e. r = 100 fm/c). Even 
the released energy jt;^'^^^'^^^'^^ drops down. This effect results from additional stopping of 
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baryon-rich fluids associated with friction with the flreball fluid. Because of this additional 
stopping, the baryon-rich fluids produce less secondary particles (and those produced are 
less energetic). Naturally, this effect is more pronounced at the energy 158 A-GeV, since 
the amount of produced secondary particles is much larger in this case than that at lower 
energies. 

At realistic values of the formation times, r = 0.5 and 1 fm/c, the effect of the interaction 
is substantially reduced. It happens because the flreball fluid starts to interact only near 
the end of the interpenetration stage. As a result, by the end of the collision process it 
looses only 10% of its available energy 

^Jreleased) ^^^^ ^ A-GcV and 30%, at Eiab = 

10.5 A-GeV. Certainly, this effect should be observable in mesonic quantities, in particular, 
in such flne observables as directed and elliptic flows. The global baryonic quantities stay 
practically unchanged at flnite r. Indeed, jr;^^^^'^^'^'^^ remains almost the same as at r = 100 
fm/c. 

The energy content of the baryon-rich fluids exceeds that of the flreball by an order 
of magnitude at Eu^ = 158 A-GeV (3440 GeV) and even more at Eu^ = 10.5 A-GeV (929 
GeV). Therefore, the interaction with the flreball fluid does not essentially change the global 
baryonic quantities. As for reflned baryonic observables, our preliminary calculations of the 
directed nucleon flow show no changes at -Eiab = 10.5 A -GeV and only slight changes in 
the mid-rapidity region at Ei^h = 158 A- GeV. This means that our previous results on the 
nucleon directed flow and its excitation function, obtained within the 2-fluid model j^, are 
not affected by the interaction between the baryon-rich and flreball fluids. 

V. CONCLUSIONS 

In this paper we have developed a 3-fluid model for simulating heavy-ion collisions in the 
range of incident energies between few to about 200 A- GeV. In addition to two baryon-rich 



fluids, which constitute the 2-fluid model 



a delayed evolution of the produced 



baryon-free (flreball) fluid is incorporated. This delay is governed by a formation time, 
during which the flreball fluid neither thermalizes nor interacts with the baryon-rich fluids. 
After the formation, it thermalizes and comes into interaction with the baryon-rich fluids. 
This interaction is estimated from elementary pion-nucleon cross-sections. 

The hydrodynamic treatment of heavy-ion collisions is an alternative to kinetic simu- 
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lations. The hydro dynamic approach has certain advantages and disadvantages. Lacking 
the microscopic feature of kinetic simulations, it overcomes their basic assumption, i.e. the 
assumption of binary colhsions, which is quite unreahstic in dense matter. It directly ad- 
dresses the nuclear EoS that is of prime interest in heavy-ion research. Furthermore, our 
3-fiuid model uses only friction forces instead of a vast body of differential cross- sect ions of 
elementary processes, which are generally unknown experimentally. Naturally, we have to 
pay for all these pleasant features of hydrodynamics: the treatment assumes that the non- 
equilibrium stage of the collision can be described by the 3-fiuid approximation. However, 
all the assumptions used are quite transparent and can be tested numerically. 

We have simulated relativistic nuclear collisions within the 3D code based on the relativis- 
tic 3-fiuid hydrodynamics combined with the EoS of the statistical mixed-phase model of 
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We performed calculations of 



the deconfinement phase transition, developed in 
nucleus-nucleus collisions Pb+Pb at -Eiab = 158 A-GeV and Au+Au at -Eiab = 10.5 A-GeV. 
To reveal the role of the interaction between fireball and baryon-rich fluids, we examined 
the evolution of global quantities of the fireball fluid. 

For zero formation time (r = 0) the interaction strongly affects the fireball fluid: It 
considerably reduces its total energy as compared to that without interaction. However, 
for realistic reasonable finite formation time, r ^ 1 fm/c, the effect of the interaction is 
substantially reduced. The fireball fluid looses only 10% of its available energy at Eig^y, = 
158 A-GeV and 30%, at i^iab = 10.5 A-GeV. Certainly, this effect should be observable in 
mesonic quantities, in particular, in such sensitive observables like directed and elliptic flows. 
Since the energy content of the baryon-rich fluids is much higher than that of the fireball 
fluid, global baryonic quantities remain insensitive to this interaction. As our preliminary 
calculations show, even the directed nucleon flow remains practically unaffected by this 
interaction. In particular, this fact justifies our previous results on the nucleon directed flow 
and its excitation function, obtained within the 2- fluid model j^. 
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